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Abstract 

The phase diagram of two bosonic chains coupled by interchain hopping and 
interchain boson boson repulsion is investigated by means of bosonization 
and renormalization group techniques. It comprises two types of charge den- 
sity waves, a conventional superfluid phase and a superfluid phase made of 
condensed boson pairs. The effect of a random potential on that system is 
also derived . Compared with the one chain case, interchain hopping strongly 
stabilizes the conventional superfluid phase with respect to Anderson localiza- 
tion. By contrast, the superfluid phase made of condensed boson pairs is much 
more unstable with respect to Anderson localization than the superfluid phase 
of the one chain system. An exact mapping of the two strongly correlated 
bosonic chains onto two coupled vortex planes gives us the phase diagram 
of the latter system without any further calculation. The vortex system ex- 
hibits two types of vortex lattice phase(one with entanglement and another 
one without entanglement) , a normal phase and a phase with bound pairs 
of vortices. The results for Anderson localization of the bosons allow us to 
study vortex pinning by columnar defects or twin boundaries and obtain the 
correlation length of the pinned vortex lattice. We show that interplane vor- 
tex hopping strongly reduces vortex pinning by columnar defects but not by 
twin boundaries. We obtain the melting temperature of the entangled solid. 
We also derive the critical current by modified Larkin-Ovchinnikov arguments 
and find a very sharp decrease of critical currents when the system moves into 
the entangled phase in the case of columnar pins. For strong columnar pins, 
there is a transition into a pinned solid phase with dislocations. We argue 
that two multicritical points may exist in the phase diagram. 

I. INTRODUCTION 

Localization effects in quantum interacting systems remain, even now an extremely chal- 
lenging problem. For fermionic systems, although interactions are necessary to describe a 
certain number of experimental situations, the non-interacting case allows for a reasonable 
starting point in most situations , allowing to deal with the complications due to disorder!. 
For bosonic systems, the theoretical situation is much less favorable since repulsive in- 
teractions have to be included from the very beginning. Despite its theoretical intricacy, 
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the problem of disordered bosons is of importance for a variety of experimental situations 
such as 4 He in aerogels or Vycori, granular superconductorsS, disordered superconducting 
films§~iand disordered Josephson junctions!!!. More recently this problem has also known a 
renewed interest^both theoretical and experimental, due to the advent of High Temperature 
superconductors". In these systems, thermal fluctuations can be strong enough to unpin 
the Abrikosov Lattice from impurities or even melt itQ0 thus restoring energy dissipation 
in the system. Finding a way to pin the vortex lattice more efficiently has therefore become 
a problem of technological importanceEl for applications of HTSC that also arises important 
questions in the basic physics of vortex systems. Columnar defects parallel to the magnetic 
field are expected to greatly enhance the pinning of vortices. From the theoretical point of 
view, it has been shown that the problem of vortex pinning on columnar defects is equivalent 
to the problem of Anderson localization in a random potential of two dimensional bosons 
with repulsive interactionsBlll. In that mapping, superfluidity of the bosonic system corre- 
sponds to a vortex liquid phase and restoration of energy dissipation whereas localization 
corresponds to pinning of vortices and absence of energy dissipation. 

The necessity of keeping the repulsive interactions in disordered bosonic systems is due to 
the fact in the non-interacting limit, the ground state corresponds to all bosons condensed 
in the lowest energy eigenstate resulting in an unphysical collapse of the system. As a 
result, one has to start with a non-disordered boson system with repulsive interactions, 
and consider the effect of randomness. In a two or three dimensional system, this is an 
extremely hard problem. One dimensional systems allow for a controlled treatment of the 
interactions therefore improving the accuracy of the subsequent treatment of disorder. It 
is thus not surprising that the first analytical solution of this problem@ was found for 
the case of one- dimensional dirty bosons. A transition between a superfluid phase and a 
localized phase upon increasing the repulsion was found. The existence of such a transition 
received numerical connrmationEl, and was extended to the more experimentally relevant 
two dimensional situation by beautiful scaling arguments^!. Since then localization of bosons 
in two dimension has been confirmed by a variety of numerical techniquesEihil. 

However despite intensive efforts a controlled microscopic theory of such a transition 
in more than one dimension is still lacking, and one has to resort to various mean field 
type approximations, the validity of which is difficult to controlii. To try to interpolate 
between the one-dimensional solution and a higher dimensional one, we thus study in this 
paper the model of a two leg bosonic ladder. Such a system has the advantage of containing 
interchain tunneling, allowing to mimic higher dimensional systems, while retaining the 
controlled treatment of the interactions possible in the one dimensional world. One thus 
expects such a study to give some clue about higher dimensional systems. For fermionic 
or spin system the study of pure@~0 or disordered§Zh§l ladders has been extremely fruitful 
in that respect. The microscopic results for the ladder can also serve as a test-bed for 
scaling arguments^! that are expected to be valid in all space dimension. Indeed, the scaling 
relations where found to be compatible with the exact exponents^ found in the single chain 
superfluid-insulator transition. In addition one expects our study to be directly relevant 
for experimental systems such as Josephson junction networks that are practical realization 
of one-dimensional bosonsSB, but often with more than one channel involved in charge 
transport. 

The ladder situation becomes also specially important when dealing with vortex systems 
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in Type II superconductor^. The superfluid to insulator transition that occurs in the dirty 
bosons system when interparticle repulsion is increased results in the equivalent vortex 
system in a transition from an unpinned vortex liquid to a pinned vortex solid as temperature 
is lowered. This pinned phase is known as "Bose Glass" phase and its physical properties 
have been analyzed in Ref. [15| and have been the object of intense experimental studiesiHl. 

Treating such system analytically proves to be difficult and possible only in the elastic 
limit where not dislocations are allowed in the lattice structured. Although it is indeed the 
case for weak disorder at low temperature, as for point like disorderly where the existence 
of a dislocation free glassy phase was proven both analyticallyil'&ii and numerically^ 
such a restriction cannot allow to approach the melting transition nor the case of strong 
disorder where dislocations obviously play a role. It is thus interesting to study special 
geometry, where dislocations could be treated exactly. This is for example the case for 
vortices confined into planes parallel to the magnetic fieldill, corresponding to coupled chains 
of one dimensional bosons. This situation may be experimentally achieved by putting 
magnetic field parallel to the CuO planes in superconductors such as YBCO or BiSSCOc 
In the case of_point like disorder this geometry allowed to confirm the existence of the Bragg 
glass phasecZH, and to study the melting and strong disorder case in a controlled way. 
However in these approaches tunneling of vortices between planes was not treated. When 
included!*! it was argued that the periodic potential induced by the layers would still lead 
to a smectic order of the vortices, with even the possibility of supersolid phases between a 
crystal at low temperature and the liquid (or smectic) at high temperatures. 

The ladder system allows for a detailed microscopic solution of such systems, taking 
both hopping, interactions and disorder into account. It is thus possible to describe more 
realistically the possible phases of a two dimensional system, and in particular to take into 
account such effects as the entanglement of the vortices. Indeed in the vortex language, boson 
tunneling translates as the ability for bosons to overcome the potential barrier that maintains 
them in the planes by gaining energy from thermal fluctuations and form "super kinks" 0. 
Also, the perturbative renormalization group arguments leading to the supersolid phase also 
applied to the two planes vortex system. Therefore, the existence of such phase can be 
tested in the two plane system going beyond perturbation theory. 

The plan of the paper is as follows: In section |J, we recall the bosonization procedure 
for bosonic systems. Then, we derive the bosonized Hamiltonian of two coupled bosonic 
chains. This allows us to obtain the phase diagram of the two chain bosonic system using 
renormalization group arguments. In section |TJ we consider the effect of random potential 
scattering. We obtain the phase diagram in the presence of disorder as well as the boson 
localization lengths and discuss the physical consequences. In section |VL we consider first the 
mapping of interacting vortices confined in planes onto interacting bosons. We give a physical 
interpretation of the quantities defined in the boson language in terms of vortices. We 
discuss the phase diagram of the vortex system in the absence of columnar disorder. Then, 
we discuss the effect of columnar disorder and compare with the one chain system. Finally, 
we give a calculation of the critical current of the coupled planes vortex system. Conclusion 
can be found in section [V|. Some technical details have been put in the appendixes. 
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II. PURE BOSON LADDER: PHASE DIAGRAM 



A. Hamiltonian 

Let us consider the most general bosonic ladder with both interchain hopping and inter- 
chain interactions. Its Hamiltonian is 



H = dx 



2m 2m 



+ J dxdyU{x - y)pBi{x)pBi{y) 



+ J dxdyU(x - y)p B2 (x)p B2 (y) + J Vp m (x)p B2 (x)dx 

+ / t±(4> Bl (x)ip B2 {x) +4) B2 (x)4) B i(x))dx (1) 



where t± is the interchain hopping, U and V are respectively the intrachain and interchain 
interactions between the bosons. Although we used for (|1|) a continuous representation it of 
course also applies to discrete models (such as e.g. the bosonic Hubbard model) provided 
the boson filling is incommensurate with the lattice (for the case of commensurate filling 
see Ref. f|6] for a single chain or Ref. |29| for a commensurate ladder). Physically, one has 
t± < and V > 0. To treat (]]]) is it convenient to use an alternate boson representation of 
the original bosonic degrees of freedom^. The salient points of this bosonization of bosons 
can be found in Appendix |^ where the single boson operator and the density operator are 
expressed as (see Appendix [A]) 

p B ( x ) = (p + ^) V e^OO-HT**) 

7T — ' 

m=— oo 

oo 

xp B (x) ~ e ld{x) \p B (x)] l/2 ~e^ } Y, e 2imWx)+7Tpox) (2) 



m=— oo 



where po is the average boson density,6'(x) = tc dx f H(x') and and n are conjugate 
operators and a is a short distance cutoff. The Hamiltonian (JJ) is bosonized by treating 
separately the two chains. In (|l|), a continuum Hamiltonian is used, but the bosonization 
procedure would also work for a lattice Hamiltonian. To do so, one introduces the two sets 
of conjugated fields $i, 4>\ for chain 1 and 8 2 , 2 for chain two. In the absence of t± and 
V the low energy properties of each chain are described by a Hamiltonian of the form ( |A1| ) 
describing the sound waves which are the phonon modes typical of a Bose superfluid. All 
the interaction effects are hidden in the so called Luttinger liquid parameters u and K (for 
each chain), u is the sound wave phase velocity whereas K controls the asymptotic decay 
of the correlation functions (see appendix [A]), u and K can be extracted directly out of 
physical quantities such as compressibility and charge stiffness, as explained in Appendix 
regardless of the precise form of the interaction U. For noninteracting bosons (U — ► 0) 
K — ► oo, whereas a local hard core repulsion would correspond to K — 1. A longer range 
interaction allows to reach all positive values of K. Taking for simplicity identical chains 
one has U\ = u 2 = u and K\ = K 2 = K. The interchain boson-boson interactions V and 
interchain hopping terms t± are expressed in terms of the fields <f>i, 6\, (f) 2 , 9 2 using (0). It is 
convenient to use the symmetric and antisymmetric linear combinations: 
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Vs = 1=—, Va = 



V2 
V2 ' 



V2 

V2 



(3) 
(4) 



In this representation symmetric degrees of freedom and antisymmetric ones decouple leading 
to the following bosonized form for the Hamiltonian ([]]) 



H = 


H a + H s 


H s = 


r dx 






U s 


1 2^ 




H a = 


r dx 




J 2~ix 











U s 



H a = u a K a (rcU a ) 2 + %-(dM 2 dxcos(V29 a ) + f dxcos{^ a ) (5) 

J 2tt L is. a 1 ira J {ira) z J 

The Luttinger liquid parameters for the symmetric and antisymmetric modes are given by 

(6) 
(7) 

where the upper sign is for the s mode. H s has the same form as ( |AT|) so that correla- 
tion functions of the symmetric degrees of freedom are easily obtained using ( |A4|) . On the 
other hand, the antisymmetric degrees of freedom are more difficult to treat. To do so one 
can exploit the analogy of H a with the Hamiltonian encountered for coupled spin chain 




system 



I. H a has two different gapped phases. When V is the most relevant field, 



acquires a mean value that minimizes the ground state energy (i.e. (<j) a ) = | for V > 0), 
whereas the correlation functions of the exponentials of 8 a all decay exponentially fast. When 
t± is the most relevant field the situation is reversed: 9 a has a mean value that minimizes 
the ground state energy (i. e. (9 a ) = 0) and all exponentials of 4> a decay exponentially. In 
the first case, the system is dominated by interchain repulsion and the density fluctuations 
between the two chains are frozen. In the second case interchain hopping dominates corre- 
sponding to the standard Josephson coupling and assuring that the superconducting phase 
is identical ((9 a ) = 0) between the two chains. The boundary K% between these two regimes 
depends on the values of V and t± . A way to estimate the boundary is to consider that the 
most relevant field is the one with the shortest correlation length. Using the standard for- 
mula for the correlation length of sine-Gordon time equations (see appendix [A]) one obtains 
for the values at the boundary (for small V, t±) 



iog(VQ 

log(*±) 



2-2K° 
2 - 1/2K< 



(8) 



For K = 1/2 the transition occurs when V = t±. Let us not that since <fi and 9 are conjugate 
fields it is not possible to get simultaneous order in both fields even if at the perturbative 
level both operators in the Hamiltonian can be relevant simultaneously. The schematic 
phase diagram is shown on Figure 0. 
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B. Correlation functions and phase diagram 



As a consequence of the existence of the gap in the antisymmetric degrees of freedom, 
most of the operators have exponentially decaying fluctuations, and the associated suscep- 
tibilities remain finite. Some operators however have algebraic correlations and diverging 
associated susceptibilities. These are the following 

a 

(9) 
(10) 

in(V20 o ) (11) 

(12) 



OsF+ 


1 6s 

= tpBi + lpB2 = e ^ 




OcDW i7rp o 


= PB1,Attpo + PB2,4Trp 




OcDW- 




e iV2(j>3 


— PBl,2irp ~ PB2,2irp 


= sin 


Obpsf 


= i>Bi{x)ipB2{x) ~ e 


i\f2Q s 



The CDW order parameters describe "crystalline" type phases, whereas the SF ones corre- 
spond to superffuid phases. The Ocdw- operator is the order parameter of an antisymmetric 
charge density wave. In that phase, there is a charge density wave on each charge of the 
ladder. These charge density waves have the same 27rp wavevector and are dephased so 
that a minimum of the density on one leg of the ladder corresponds to a maximum of the 
density on the other leg. Since we consider mainly interchain repulsion V > we confine 
ourselves to this order parameter. For an attractive interchain interaction one would have a 
corresponding Ocdw+ order parameter (with cos(v^20 a ) instead of sin(y / 20 a )) correspond- 
ing to a phase for which the two 27rp CDW on each chain would be in phase. The O CDW ^ P0 
operator is the order parameter of a symmetric charge density wave with 4-7rpo wavevector. 
This corresponds to bosons that form a charge density wave with a period half the one of 
the CDW-. The physical picture is that of bosons forming a CDW- but hopping from one 
chain to the other thus minimizing the hopping energy. This has the effect of halving the 
period and changing an antisymmetric charge density wave into a symmetric one. Strictly 
speaking, the order parameter for the 47rpo charge density wave is in fact: 

Ocdw^o = e lVE ^ cos V8<f) a (13) 

and expression ([II]) results from the combination of the term V cos \fS<p a of equation @ 
with equation ( |T3"D in the perturbative expansion of the correlation functions. The Osf+ 
operator is the order parameter of a superfluid phase in which there is Bose condensation of 
the bosons that occupy the bonding band. This is the normal superfluid order parameter. 
On the other hand the Obpsf operator is the order parameter of a superfluid phase in 
which pairs of bosons that are formed along the rungs present a Bose condensation. In a 
lattice system, these pairs would be formed of a boson on i-th site of chain 1, and another 
boson on the (i + l)-th site of chain 2 for a repulsive V (on the same site for V < 0). 
This implies that the BPSF also contains antisymmetric charge density wave (i. e. CDW-) 
fluctuations. Which phase is realized depends on whether interchain hopping or interchain 
interactions dominate. The CDW 47rpo and SF+ phases are typical of a system dominated by 
interchain hopping whereas the CDW- and BPSF phase correspond of a system dominated 
by interchain repulsion. Using (^) one can obtain the generic phase diagram for the boson 
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ladder. The gapped fields, their mean value and the phase that obtains are listed in table |. 
The resulting phase diagram is shown on Figure |2|. For moderate intrachain repulsion 
K s > 1, and only the superfluid phase is possible. If V is small with respect to u/K, 
K a = K + 0(KV/u), K s = K + OiKVju) and one moves on a special line of the phase 
diagram with K a ~ K s (dotted line on Figure 0). Note that, as can be seen from (||), 
increasing the small q interchain repulsion pushes paradoxically the system to a regime 
where the interchain hopping is more relevant and thus destabilizes the crystal. This is due 
to the fact that now regions of high density on one chain face regions of low density on 
the other making it easier for the hopping to act. Therefore, in a system of two coupled 
chains of bosons with repulsive local interchain interactions the BPSF does not exist (see 
Figure |), which is physically reasonable. Instead, one has a superfluid SF+ phase for 
K > 1/2 and a CDW- phase for K < 1/2 corresponding to a boson "crystal". For attractive 
enough interchain interactions one could expect to stabilize the BPSF phase. Finally let 
us make the connection between the bosonic ladder and other models. First using the 
standard equivalence between spins and hard core bosons (see e.g.E§), the bosonic ladder 
can be mapped to a spin ladder. In the absence of magnetic field for the spin system, 
one would have a commensurate filling of one boson every two sites. The massive phases 
(Haldane or singlet) for the spin system would thus correspond to an insulating Mott phases 
for the bosonic ladder, where bosons are frozen along the chains but able to hop between 
the two chains. The Ising antiferromagnetic phase corresponds to the boson "crystal" . The 
massless XY phases are the superconducting ones, the BPSF phase being equivalent to the 
XY2 phase. Since such systems have been studied in details both for the pure system and 
in the presence of disorder we refer the reader to Ref. ^ for a more detailed analysis of 
the physical properties. Here we mainly focus on the incommensurate situation that can 
be mapped to the spin ladder in the presence of a finite magnetic field, playing the role for 
the bosonic ladder of the chemical potential. Such bosonic ladders can also be mapped to 
the standard fermionic ladder in the presence of an on-site attractionil. In that case the 
electrons pair in bound states and single electron hopping becomes irrelevant. The system 
becomes equivalent to pairs of carriers hopping for site to site with an effective hopping of 
the order of t 2 /\U\ for a strong attraction U. The interchain hopping can also be viewed in 
that the standard Josephson coupling between two superconducting chains. More 

details can be found in ref. |28|. 



III. EFFECT OF DISORDER 

Let us now add disorder to the bosonic ladder. In a continuum system, the only possible 
form of disorder is a random chemical potential. In the case of a lattice system, another 
possible form of disorder besides a random on site potential is random hopping along the 
chains. Although these two types of disorder can lead in principle to different physical phases 
when particle hole symmetry exists (as is for example the case for spin chainsH^H), for an 
incommensurate density, they lead to the same coupling in the continuum limit. Since we 
confine our study of the bosonic ladder to such a case, we can restrict ourselves to consider 
a random chemical potential. The disorder is thus modeled as two uncorrelated random 
Gaussian potentials in chains 1 and 2. In the continuum limit, the coupling of the bosons 



7 



with those random potentials is: 

-^disorder = J dx [Vi(x)p B l(x) + V 2 (x)p B 2 (x)] (14) 



where V p (x)V p i(x') = D5 P:P '5(x — x'). Upon bosonization the full Hamiltonian including the 
coupling with disorder has the following form 

H = H a + H s + -^disorder 
^disorder = fdxJ2 [V?™ P ° ( x ) e * m (v / 2(^+<MM) + yTmxpo ^ym^^-^x)) + R ^ ^ 



where H a ,H s are given by (0) and V P 2rnvpo (x)V pl mnpo (x') = D m 8{x — x')8 P)P i. The q ~ part of 
the random potentials has been dropped since it does not affects the gaps for small disorder 
and does not induce localization as discussed in Ref. As usual the most relevant coupling 



to disorder is given by the lowest harmonic0 so we retain for the moment only V 27Tf>0 and 
D± ~ D. Higher harmonics can be taken into account in a similar way using (0L The 
analysis of the phase diagram of (|i~5|) proceeds in a similar way than for spin ladderso, with 
the important physical difference of a gapless symmetric sector for the boson case. Disorder 
has different effects on the phases dominated by interchain hopping (i.e. with (8 a ) = 0) 
corresponding to sectors I and II of table | and the ones dominated by interchain repulsion 
(i.e. with (cf) a ) = corresponding to sectors II and IV of table |I[ In each case, two phases 
are in competition. In the former case, there is a competition between a 47rpo CDW and 
a SF+ phase. The 4irpo CDW phase can gain energy by pinning on the 4irpo component 
of the random potential whereas in the SF+ phase, density fluctuations are too small to 
permit gaining energy from the random potential. In the latter case, the competition is 
between the 2irpo charge density wave and the BPSF phase. The mechanism is similar, 
however one expects that the stabilization of the CDW- is stronger than the stabilization of 
the 47rpo CDW. Technically, the case of the competition of the CDW- and the BPSF phases 
is simpler and is discussed first. 



A. System dominated by interchain repulsion 

When the system is dominated by interchain repulsion, the two phases in competition 
are the CDW— and the BPSF phase. The effective Hamiltonian describing the coupling 
with disorder is obtained by replacing (f> a by its mean value in fllpp. It has the form : 

^effective = / V eS .(x)e^ (16) 

where V e s,(x) =\ (e lv ^^ a ) \ i(Vi(x) — V 2 {x)) and only the most relevant (in the RG sense) 
terms have been kept. From (|1^) the renormalization group equation for disorder is straight- 
forwardly obtained in the form: 

^ = &-K.)D (17) 
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The renormalization group equation for K is : 

dK „ 9 

where C is a non-universal prefactor. In a similar way than for a single chain0 these RG 
equations have a fixed line on which disorder is irrelevant, and a strong coupling disordered 
fixed point. For small D, the renormalization of K can be neglected and (|17D is sufficient 
to determine whether the system flows to the fixed line or to the strong coupling fixed 
point. According to ([17])) the disorder is relevant for K s < 3 and irrelevant for K s > 3. 
The strong coupling fixed point is in this case an antisymmetric charge density wave CDW- 
pinned on the 2npo component of disorder@@. Therefore, the BPSF phase is unstable with 
respect to a very small amount of disorder unless K s > 3. This is due to the coupling of 
the random potential to the slowly decaying CDW- fluctuations (see equation ([□])) that are 
subdominant in the BPSF phase . When interchain density- density interactions dominate 
the superfluid-insulator transition occurs thus for much less repulsive interactions (K = 3) 
than for the single chain case (K = 3/2). The situation is reminiscent of what happens in 
fermionic ladders^, and is due to the reduction of quantum fluctuations by the freezing of 
the transverse degrees of freedom. The localization length of the bosons in the pinned phase 
is easily extracted from the RG equation as the length for which the renormalized disorder 
is of order one. It gives 



/ 1 \ 3-K s 

iloc.,2ch. ~ (jj) (19) 

and as can be expected is much shorter than the localization length of a single chain 

/ 1 \ 3-2K s 

^loc.,lch. ~ ( J (20) 

Note that as for the fermion caseH the change in localization length occurs through a change 
in exponent and not as a simple multiplicative factor. Since the Bose fluid that condenses is 
made of bosonic pairs, to obtain the exponent of the correlation function of the superfluid 
order parameter at the transition one must consider 

(T T O j BPSF (x,r)O BPSF (0,0)) = ( 1 V ' (21) 

\X l + (UT) 2 J 

as the superfluid order parameter correlation function. That correlation function decays as 
r -i/K s j n superfluid phase and at the transition point, K s = 3. Thus, at the transition 
point T) c — 1/3. It is interesting to remark that the same exponent at the transition point 
obtains in the case of one bosonic chainS. Using the RG one can also compute the conduc- 
tivity in the case of a charged Bose fluid made for instance of Cooper pairs. This gives an 
indication of the amount of dissipation that is generated in the system upon introduction of 
disorder. In a neutral boson system, conductivity measures the normal current that would 
be induced by a gradient in chemical potential. The high frequency dependence of the a.c. 
conductivity can be obtained in the form : 
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(22) 



From the a. c. conductivity, it is easy to deduce 
c. conductivity. In the calculation of Refs. [L6IJ2T 



using the methods exposed in Refs. |16|,P 
the temperature dependence of the d. 
frequency just serves in determining a cut-off length — . Similarly, temperature provides a 
cutoff length This leads to the following expression for the temperature dependence of 
the d. c. conductivity : 



a d . c .(T)=T 



2-K a 



(23) 



For T <C T pin . = p 2 - or uj <C T pin ., the conductivity goes to zero exponentially fast. 



B. System dominated by interchain hopping 

In a system dominated by interchain hopping, the two phases in competition are the 
SF+ phase and the CDW 47rpo . Since in a system dominated by interchain hopping, (6 a ) = 0, 
the massive degrees of freedom in ([15]) must be integrated out in order to have a non trivial 
effective coupling to disorder. Alternatively one can couple to a higher harmonic of the 
random potential. Using the method of Ref. [28|,the effective coupling to disorder gives 

#imp, off. = / dxV eS .(x)e^ + H. c. (24) 

where V P & (x)V P fe (x') = CiD 2 5(x — x'). C\ is a constant that comes from the integration of 
the gapped modest . Physically, this represents the coupling of the impurity potential with 
the 47rp component of the density (see (0)). Note that the genuine 47rp component of the 
random potential is coupled to a cos(a/80 q ), so that its contribution is a second order term 
V 2 D 2 . In this section we will thus use D for the effective disorder. The RG equation for 
disorder is obtained from fl2"4"|) as 

^ = (3 - 4K S )D (25) 

The impurity potential is relevant for K s < | and irrelevant for K s > | . Therefore, the 
SF+ phase is stable with respect to Anderson localization for K s > |, corresponding to 
extremely repulsive interactions (this would need hard core bosons plus a nearest neighbor 
repulsion). For K s < |, the randomness is relevant and the SF+ phase is replaced by a 47rp 
CDW pinned on the disorder. The pinning length of the CDW, which is also the localization 
length of the bosons, is given by : 

(1 \ 3-4if s 

This time, the localization length is much larger than the localization length of the one 
chain system (pD"D, indicating that the two chain system is less localized than its one chain 
counterpart. In particular, a system of two hard-core bosons chains (which has K = 1) with 
interchain hopping would be superfluid even in the presence of a weak random potential, 
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whereas its one chain counterpart would have been already localized (K c = 3/2). The 
critical exponent of the superfluid order parameter correlation function is obtained by the 
method of the preceding section. In the superfluid phase, the 2 point correlation function 
of SF+ decay as r ~ l / 4Ka . Since the phase boundary is at K s = 3/4, the exponent of the 
superfluid order parameter at the transition point is 1/3 as in the preceding section. The 
high frequency dependence of the conductivity is computed as: 



0" a 



CO 



4:K S ~A 



(27) 



using the methods of Refs. [Iq , |28| . Due to the equivalence of frequency and temperature one 
also has: 



cx d . c .(T) = T 



2-AK s 



(28) 



As in the preceding section, high frequency conductivity refers to a.c. conductivity at 
frequency above the characteristic pinning frequency w,s//2ch.- The full phase diagram in the 
presence of disorder is represented on Fig. |3| 



C. Physical discussion 

Although most of the physical realization of bosonic ladders are likely to come from 
classical systems (see next section), its study also present direct interest for quantum sys- 
tems. Indeed recently practical realization of one dimensional bosonic systems have been 
obtained, mostly based on Josephson junction arraysi. It would be interesting to check in 
these systems various of the predictions made here. The easiest to determine is certainly 
the strong stabilization of the superfluid phase for the ladder in comparison with the single 
chain case. For the pure system, for a commensurability of one boson per site the Mott 
transition occurs at K — 1 (the hard core limit) for the ladder against K = 200 for a 
single chain, i.e. much less repulsive interactions. In the presence of disorder a similar sta- 
bilization occurs since the transition goes from K = 3/2 for a single chain to K = 3/4 for 
the ladder. Maybe even easier to investigate would be the bosonic ladder in the_presence 
of a single barrier. The properties of uniform disorder and a single impurityHH can be 
easily connected^ which allows from the present study to get the parameters for which the 
barrier would remain transparent. For a single bosonic chain the barrier is not effective for 
K > 1 whereas for the ladder the limit would be K > 1/2. Another, albeit more theoretical 
interest of the bosonic ladder is to provides another analytic system on which the general 
scaling arguments^! proposed for the superfluid-insulator transition can be checked. Indeed 
at that time these arguments could only be benchmarked on the only existing analytical so- 
lution for a superfluid-disordered transition, namely the d — 1 single chain cas Although 
since then thesescaling predictions have been found in agreement with various numerical 
simulations00@, it is still interesting to verify them on models for which the exponents can 
be computed exactly. The equality z = d for the dynamical exponent clearly still holds for 
the ladder since the system is still described by a Luttinger Liquid. Similarly it is easy to 
check from the RG equations that the compressibility still remains finite at the transition, 
again in agreement with the scaling prediction. Finallyin Ref. [TBI, a bound for the exponent 
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of the superfluid order parameter at the transition point has been derived by combining 
an exact inequality with scaling. When specifying to z = d = 1, that bound is r\ c < 1. 
Since r\ c = 1/3 at the transition, both in the case of the BPSF to insulator transition and 
SF+ to insulator transition, the bound is indeed verified for the ladder. It would clearly be 
interesting to extend such an analysis to ladders with more than two legs to complete the 
check. Finally the two leg ladder prompts for the interesting question of the existence of a 
universal exponent for the superfluid transition in one dimension rj = 1/3, regardless of the 
number of legs. Although a detailed check of this statement would need a complete solution 
of the models with more than two legs, simple arguments given in Appendix [FJ suggest that 
this is indeed the case at least in systems in which all modes but the symmetric one are 
gapped. This is the case of all systems with an even number of legs. However, this lets the 
possibility that the exponent could be different from 1/3 in systems with an odd number of 
chains. An interesting question is whether the exponent is also 1/3 for the transition from 
the pinned charge density wave that occurs along the line 3/4 < K s < 3, K a — 1/2 or if it 
depends of K s . It is difficult to answer this question in the framework of the approximation 
of a disorder much smaller than the gaps, since such approximation breaks down precisely 
in the vicinity of that line. 

Let us note that of course the universality of the exponent holds here for the case where 
the disorder is much weaker than all the gaps in the system. The argument cannot be 
applied directly for a two dimensional system. 

IV. COUPLED VORTEX PLANES PROBLEM 

Although specially engineered quantum systems are necessary to have coupled bosonic 
chains, such systems have a natural realization in the classical world of vortices. Indeed, 
using a standard mappingBEl space time trajectories of bosonic particles can be mapped 
onto a problem of interacting flux lines. The disorder in the quantum problem is time (i.e. 
z) independent , corresponding in the vortex problem to columnar defects. Defects of that 
type can be artificially produced by irradiation. In the presence of such defects the vortex 
system is known to give rise to a "Bose glass" phase corresponding to the localized phase of 
two dimensional bosons. On the other hand a very interesting situation occurs in a layered 
system if the magnetic field is aligned along the planes direction. In that case it has been 
proposed^ that the additional periodic potential due to the planes can confine the vortices 
and lead to a smectic like state. If the confinement is strong enough one can view the vortices 
as existing only in planes with repulsive interactions between adjacent planes. Such coupled 
planes models are especially interesting to study in their own since they allow a detailed 
description of the melting transition for thepure system as well as a study of the relevance 
of dislocations when disorder is presentS™. However, in most of the studies of such 
systems hopping of the vortices between planes was neglected. When hopping is included, 
it was proposed^ for the pure system that three phases would exist: a solid phase where no 
hopping is present and the system crystallize due to interplane repulsion, an intermediate 
supersolid phase where both the interchain repulsion and the hopping are relevant, and a 
high temperature liquid phase. Since the model does not allow for vortices to exist between 
the planes the liquid phase also possesses a smectic order along the z direction. For an 



12 



infinite number of planes such phases can only be determined by looking at the relevance of 
the various operators (interchain repulsion and interchain hopping). For two coupled vortex 
planes we can use the results of the previous sections to give precise answers to the nature 
of the various phases. Since the perturbative study of two planes and an infinite number of 
planes are identical we expect the results to be also of relevance for the more general case. 
In addition one can study the effect of columnar defects of twin boundary disorder on the 
phases of the pure system. The more delicate case of point like disorder will be considered 
elsewhere. The particular case of two coupled planes of vortices could also be of relevance to 
artificial layered structures where superconducting structures can be separated by arbitrary 
heights of insulating material&@. 



A. Mapping of the interacting vortices problem onto a bosonic problem 

1. Path integral reformulation of the interacting vortices problem in terms of interacting bosons 

In this section, we consider a bilayered superconductor in a magnetic field H > H C1 
parallel to the bilayers. We make the approximation that the bilayers are sufficiently far 
apart so that the Josephson coupling of the bilayers can be neglected. We suppose that the 
bilayers exert a strong enough attraction on the vortices so that they remain confined in the 
planes. We also suppose that there are columnar defects parallel to the planes and parallel 
to the magnetic field as shown on Fig. Such system may be realized with copper oxide 
superconductors with coupled planes (such as BiSSCO) or by epitaxial growth of artificial 
structures&~0 with layers of two superconducting materials, one of them being more efficient 
in confining vortices. 

In the following, we design by Oz the direction of the columns and of the magnetic field, 
by Oy the direction perpendicular to the plane and Ox the direction such that (Ox, Oy, Oz) 
forms an orthonormal basis. The position of a point on a flux line i of coordinate z along 
the Oz axis is given by r,(z) = (ri x (z) , ri y (z)) . The Hamiltonian for that system is then: 

H = fdz T,(P) (z) + lj2U(v i (z)-r j (z))+Y:(V(r iy )+V(r iy -2a))+j:v(v i (z)) 

i \ / iy^j i i 

(29) 



(3 being the height of the planes, 2a their distance, U being the vortex-vortex repulsive 
potential, V the vortex plane interaction and v being the random potential induced by the 
columnar defects. We consider v to be small, corresponding to weak columnar disorder 
parallel to the applied magnetic field. Although experimentally such a situation is more 
difficult to achieve than strong pinning, it allows for an analytic treatment of the problem 
and at least give a qualitative understanding of collective pinning in such systems. This 
collective pinning may be probed in high field experiments. 
The partition function of that system 

VTi(z)exp(~) (30) 
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is identicalcilO to the one of a system of interacting bosons at temperature j3 with a 
Planck's constant T. Assuming that the attractive potential imposed on the vortices by 
the planes is strong enough, the tight binding approximation can be used to simplify the 
boson Hamiltoniaril3. The calculations can be found in appendix 0. The resulting Hamilto- 
nian is identical to (fl|), i.e. two bosonic chains, with interchain and intrachain boson-boson 
interactions, interchain hopping and a random intrachain potential. Therefore, we only have 
to translate the results already derived for the bosonic system in section |11 A| in the vortex 
language. 

2. Interpretation of the bosonized vortex Hamiltonian in terms of vortices 

In the case of a single chain, the field <f>(x) has a simple interpretation in terms of flux 
lines. If u(x) represent the displacements of the vortex lattice in the elastic theory , one 
has: 4>(x) = —irp u(x) by identifying the expressions for the q ~ component of the 
density. As a consequence, the action of the bosonized theory identifies with the elastic 
Hamiltonian of the vortex lattice. The compressibility x an d the charge stiffness T> of the 
bosonic system defined respectively in equations Eqs. (|A5|) and ([AT]) have also a simple 



physical interpretation in terms of vortices. x~ l is proportional to the bulk modulus of the 
vortex system. x _1 is also proportional to the magnetic susceptibility of the system since the 
vortex number is proportional to the magnetic field. 1/T> is proportional the shear modulus 
of the vortex lattice. This results from writing the action of the bosonized theory (which 
is also the Hamiltonian for the vortex lattice elastic theory) and considering the effect of 
making d z <f) — > d z <p + A on the energy. A detailed derivation can be found in appendix [E]. In 
the vortex language, the expression ( |A3| ) has been derived directly in Ref. |35] starting from 
the elastic theory of the one dimensional vortex lattice. However, in the case of a two plane 
system, the mapping on bosons simplifies the description of vortex tunneling. 

B. Phase diagram of the vortex system 

1. Flux lines in the absence of disorder: 

In the absence of disorder, the vortex system is equivalent to the bosonic system in the 
absence of random potential discussed in section |il A]. The parameter K s controlling the 



stability of the different phases depends on temperature, vortex density and elastic constants 
of type vortex lattice as shown in Appendix [E| Its expression is: 

(Cud) 1 / 2 

Each of the four bosonic phases translates into a phase of the vortex system. In terms 
of flux lines, the CDW- is a crystal phase with a periodic array of flux lines in each plane. 
The periodic array of vortices in the second plane is translated of half a period along the x 
axis with respect to the one that lies in the first plane in order to minimize the repulsion 
energy. In that phase, there is no interplane hopping of vortices. The CDW- is a strict Id 
analog of a vortex lattice. 
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By contrast, the CDW 47rpo has the same disposition of vortices but they are constantly 
hopping from one plane to the other leading to an effective period of the vortex lattice that 
is half the one of the CDW-. This is a vortex lattice with "super kinks" from one plane to the 
other. Due to the hopping between the planes this phase is the equivalent in this restricted 
geometry of an entangled solid phase. 

The SF+ is a vortex liquid phase in which a flux line is making "superkinks" from one 
plane to the other. In this phase, vortices are unpinned. This implies strong dissipation in 
the original vortex problem (i.e. in the vortex system superconductivity is lost). Since in 
this phase hopping is relevant vortex entanglement is expected. In other words, the SF+ 
phase maps onto an entangled liquid. 

On the other hand, the BPSF phase is a vortex liquid in which flux lines that are in 
different planes are bound to each other. No hopping between planes takes place. 

As usual in a two dimensional classical system, these phases do not have true long range 
order but algebraically divergent correlations. It is generally expected that by coupling 
weakly a lot of two chains systems, the phase that would exist would be the one with the 
most divergent fluctuations in the isolated two chain system as suggested by mean field 
theory. In the case of a coupled chain system, there would be true superfluidity of the boson 
system (leading to an exponential decay of density fluctuations i.e. a true vortex liquid) or 
to formation of a vortex lattice with true long range order. 

In the case: V(y) = — V 5(y) and U(r) = U 5(r), t± is much larger than V, so that 
K s ~ K a ~ K . This case corresponds to the dotted line on Fig. |2|. 

As a consequence, this system has 2 phases: a CDW- for K s < 1/2, and the SF+ phase 
for K m > 1/2. We assume here V ~ t± for simplicity otherwise the boundary is displaced 
according to Figure [l|. The behavior of vortices in these two phases is represented on figure 
||. In terms of vortices, there is a critical temperature (see appendix ||) T rn = — such 
that for T < T m the vortex system has strong fluctuations towards the vortex lattice state, 
and for T > T m strong fluctuations towards the vortex liquid state as shown on figure |6|. 

When vortex density (i.e. applied magnetic field) is increased, the "melting temperature" 
T m decreases. Also, the larger the elastic constants, the higher the melting temperature. 
Such results are very similar to the higher dimensional phase diagram. 

Note that since the and 9 fields are conjugated an intermediate "supersolid" phase 
where one would have both a crystalline and a superfluid order, cannot exist in this system, 
even if perturbatively both operators would be relevant for 1/4 < K < 1. A more detailed 
proof is given for the special temperature corresponding to K = 1/2 in Appendix [F[ Whether 
such supersolid phase can exist for an infinite number of planes as advocated in ref. pE5] 
remains an open question. 

2. Effect of columnar disorder on the vortex system 

In the presence of columnar disorder, the vortex system at finite temperature is equivalent 
to the two chain boson system in the presence of a random potential discussed in section 

If the interplane interactions are not too strong K s ~ K and the various phases can be 
read from figure The pinning of vortices is then controlled by the parameter K s = 
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where T m is the melting temperature. 

The main effect is to replace the solid-liquid transition that existed in the pure system at 
T = T m by a double transition. The solid phase is replaced by a pined solid corresponding to 
the 27rp vortex lattice phase. This phase is a Bose glass. At the pure melting temperature 
T m , the solid phase disappears and interplane interactions become irrelevant. The Bose glass 
phase is replaced by a pinned CDW 47rpo vortex glass phase. Since this phase has dominant 
symmetric density fluctuations is possesses some (disordered) solid type order. On the other 
hand interplane hopping being relevant, this phase could also be viewed as a pinned solid 
with entanglement, with intermediate properties between a solid and a liquid. It is only 
above a higher temperature T L = 3/2T m that this entangled solid disappear to give place to 
the genuine liquid, where disorder is irrelevant, as shown on figure [7|. As is expected in the 
presence of columnar disorderll the melting temperature is renormalized upward. 

In the two plane vortex system, the pinning of the 2irp vortex lattice is much stronger 
than the pinning of the 47rp vortex lattice, as can be shown by the calculation of the pinning 
lengths and critical currents detailed in the following section. The reason is that the 47rp 
vortex lattice is destabilized by formation of superkinks that permit to average the disorder 
and reduce its efficiency in pinning the vortex system. This superkink formation is possible 
only when temperature is large enough to overcome the deformation energy that is needed 
to bend the vortices, i.e. T > T m . In the 27rp vortex lattice on the other hand, there is an 
enhancement of pinning due to strong vortex repulsion that makes the vortex system much 
more sensitive to the antisymmetric component of the random potential. For finite disorder 
since the pinning of the 2wpo vortex lattice is much stronger than the one of the 47rpo vortex 
lattice the transition line will be shifted upward from K a = 1/2, i.e. to temperatures higher 
than T m . Columnar defects obviously tend to make hopping between planes less relevant. 
However, we are in the limit where hopping dominates over columnar disorder, so that the 
difference in columnar pinning potential between the planes is simply averaged out. Such 
effect appears nevertheless in renormalization group treatment around the decoupled chains 
fixed point. 

It is interesting to note that if the disorder is perfectly correlated between the planes, as 
for example would be the case for twin boundaries, it can only couple to the 47rp component 
of the density in both the 27rpo vortex lattice and the 47rp vortex lattice, leading to the same 
pinning lengths in the two phases up to prefactors. This is easily seen by putting V\ = V2 in 
Eq. (|15D . Such result is due to the fact that the disorder being the same in the two planes, 
it does not couple to the antisymmetric density of the vortex solid phase, and the averaging 
effect of vortex hopping cannot lead any more to a reduction of vortex pinning. 

As we have seen above, the smaller the stronger the pinning. Thus, at a given 
temperature, a stiff vortex lattice having a higher T m is much more easily pinned than a 
softer one. Also, by reducing vortex density (i.e. reducing the applied magnetic field) one 
increases T m and thus the pinning of the vortex lattice. 

C. Physical discussion 
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1. Weak disorder limit 



The two planes vortex system shows two superconducting phases with different vortex 
pinning lengths. The first one, the 27rp vortex solid, is dominated by interplane vortex 

repulsion and has the shortest vortex pinning length : l 27Tpo ~ 3 ~ As ( see Eq. (pOD). 

Expressed in terms of temperature, this gives h-npa ~ (t)) 3<Ti T \ where T\ = QT m . It is 
a Bose Glass phase, and at the melting point, the vortex pinning length becomes infinite 
with an essential singularity. Such singularity is typical of a two dimensional system. The 
other superconducting phase, the 2npQ phase is dominated by formation of superkinks and 

has a larger pinning length given by Eq. (|26|) Z 47rpo ~ 3 ~ iKs since the pinning potential 
is averaged out on the planes. In terms of temperature, the behavior of the pinning length 

in this second phase is Ia-k Po ~ (5) 3(Tl_t) , where T L = §T m . It is again a Bose Glass phase, 
but the underlying lattice has half the period of the underlying lattice of the previous phase. 
Similarly to the previous phase, the correlation length shows an essential singularity at the 
melting point. These two phases could be probed in neutron scattering experiments, and 
could be distinguished by measuring the wavevector of the underlying lattice. 

A quantity of interest for vortices is of course the critical current. A proper calculation 
of the critical current involves a full solution of the dynamics of vortices whereas all our 
calculations are restricted to thermal equilibrium. Nevertheless, a simple estimate of the 
critical current along the lines of the Larkin Ovchinnikov argument can be made using 
only the equilibrium results. In the Larkin-Ovchinnikov theory, R c designates the distance 
beyond which the random force approximation breaks down. It is argued that for r>i? c the 
vortices are pinned as independent bundles. In the Larkin-Ovchinnikov theory, the balance 
between random force due to disorder and the force F c due to the critical current implies 
that : 

(32) 

for vortices confined to planes parallel to the magnetic field in the presence of columnar 
disorder. 

The distance beyond which vortices are uncorrelated is just the localization length of 
the bosons. Replacing R c with the boson localization length gives a power law dependence 
of the critical force on disorder. It can be checked first that in the limit T — > one 
recovers the Larkin-Ovchinnikov result and second that the critical current goes to zero at 
the localization-delocalization transition for the bosons indicating the onset of dissipation in 
the vortex problem. In the pinned 27rp vortex solid phase the critical current is given by: 

8T m -T 

F c ~ £)2(6T m -T) (33) 

In the pinned 47rpo entangled vortex solid phase one has: 

2T m -T 

F C ~D' AT ™- 2T (34) 
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The 4irpo Bose Glass does not extend down to OK, so that there is no reason to expect 
recovering the Larkin-Ovchinnikov result for T — > 0. For T < T m it is replaced by the 27rp 
Bose Glass as explained in section |1V C 3| . At T = T m , there is a transition from the 27rp 
vortex liquid to the 47rp vortex liquid. At that point (see figure |8|), the critical current drops 
from F c ~ .D 7//10 for T — > to F c ~ D for T — > T+. This apparent discontinuity is a 
consequence of the approximations made. However, one expects a very rapid drop of F c as 
T traverses T m . 



2. Qualitative discussion of the strongly disordered case 

In the strongly disordered case, disorder can localize bosons before interchain hopping 
or interchain repulsion can open a gap in the antisymmetric degrees of freedom. In the case 
of vortices, this corresponds to independent pinning of vortices in each plane with neither 
interplane repulsion of vortices nor hopping of vortices. In such case, the two planes are 
completely decoupled and there is no correlation of vortex positions between the two planes. 
This phase can be seen as a pinned solid containing dislocations between the planes. This 
transition is the analogous of the transition between a Bragg glass phase and a vortex glass 
one occurring for point like disorderEE. In the case of correlated disorder it corresponds to 
the transition between a dislocation free Bose glass (or Bragg Bose glassll) to the "standard" 
Bose glass containing dislocations. The condition for existence of that phase is that the 
pinning length of vortices in one plane is much shorter than both the correlation length due 
to interplane repulsion and the correlation length due to interplane hopping. This leads to: 

3T m -T 

D > V^=t (35) 

T(3T m -T) 

D^>t ± 2 t - t ™ (36) 
The resulting phase diagram is drawn on figure |9|. The transition line between the phase with 

3T m -T T(3T m -T) 

"superkinks" and the Bose Glass without dislocations corresponds to V 2T ™- T = t± 2T ~ T ™ . 

T(3T m -T) 3T m -T 

This vertical line ends at a point D = t± 2T - T ™ = V 2T ™- T . This suggests a possible mul- 
ticritical point separating a Bose glass phase with dislocations, a Bose glass phase without 
dislocations (i.e. a Bragg Bose Glass) and a Bose Glass phase with "superkinks" . 

For T > 3T m , the Bose glass phase with dislocation disappears and is replaced by a liquid 
phase. The transition line between the pinned entangled solid and the liquid phase starts 
at T = |T m for D — > 0. The melting temperature of the pinned entangled solid increases 
with disorder since the pinned solid gains energy from the random potential with respect to 
the liquid. The melting line of the Bose glass with dislocations can cut the melting line of 
the Bose glass with "superkinks" . The intersection may correspond to a second multicritical 
point. The two possible multicritical points lie in the shaded zones of Fig. ^. It is possible 
to obtain the critical currents by the method of the previous section. One obtains: 

4T m — T 

F c = D^^t (37) 
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3. comparison of the one plane and the two plane vortex system 



For the two planes vortex system, pinning depends on two factors. The first one is of 
course temperature. At high enough temperatures, vortices always depin. The second factor 
is the relative magnitude of t± and V. t± depends on the potential that confines vortices 
in the planes but also on the temperature via the unperturbed wavefunction (see appendix 
for the case of a 5 potential). V depends both on temperature and vortex- vortex repulsion. 
It also depends on the potential that confines vortices in the plane via the unperturbed 
wavefunction. Therefore, the phase diagram of the non-disordered two planes system as 
a function of temperature is strongly affected by the potential that are chosen to model 



vortex-vortex interactions and vortex planes interactions. It is known from section [IV B 1 
that in the absence of disorder the two plane vortex system is in a vortex lattice state for 
T < T m and in a vortex liquid state for T > T m . Also, from section |I V B 1| we know that 
for T < T m , we have a 2np Bose Glass, for T m < T < T L , a pinned 47rp vortex lattice 
containing superkinks, and for T > Tl vortices get unpinned. 

Remarkably, the two temperature are related by the simple relation T L = |T m . It is also 
interesting to remark that at T = T m there is a discontinuity in vortex pinning length and in 
critical current. Of course, such a discontinuity is the result of the neglect of disorder with 
respect to the gap, and a more exact treatment would lead to a continuous albeit extremely 
rapid decrease of the critical current with temperature. 

For the one plane system on the other hand, there is only one Bose glass phase and 
vortex depinning is solely driven by temperature. The transition Bose Glass to vortex liquid 
occurs at a higher temperature 7\ p i ane = 2T L and for given temperature and given amount 
of disorder the vortex pinning length is shorter in the one plane system if T > T m . On the 
other hand, for T < T m , the pinning length of the two planes system is the shortest due 
to the transition from the 4irpo to the 27rpo Bose glass as the temperature is lowered. As 
a consequence, for T < T m vortex pinning is the strongest in two planes systems whereas 
for T > T m it is the strongest in one plane systems. What distinguishes sharply one plane 
vortex systems from two planes vortex systems is the existence of a two step transition from 
the Bose glass phase to the vortex liquid. The first step is the formation of superkinks that 
halve the periodicity of the underlying lattice, leading to a second Bose glass phase much 
less well pinned by columnar disorder, and in a second step the melting of the new Bose 
glass phase. 

All these considerations are valid for T > 0. However, for T = 0, thermal fluctuations 
are not available to screen the pinning potential. It is therefore interesting when comparing 
the one plane and the two plane system to consider thelimit T — > 0. In the T — > limit, it 
is possible to obtain the pinning length by an Imry-MaE3 type argument. One requires that 
the phase of the field coupled to the random potential jumps by 2ir on a length L. This 
leads to a variation of <p s on the length L of 6(f) s ii^/2. The resulting pinning energy is: 

E(D = 2j ^-\l^ (38) 
Minimizing with respect to L gives: 
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-^pin.,2ch. 

For one plane, the same procedure gives : 

Lp'm., 1 ch. 



167^c|\ 1/3 



(39) 



(40) 



The pinning length is therefore larger in the two plane system by a factor 2 1 / 3 than in the 
one chain system. 

Therefore, at T=0, we expect two regimes: For Z) 1 / 3 3> V 1 ^ 2 , the two planes are mutually 
decoupled and individually pinned on disorder. For D 1 ^ 3 <C V 1 ^ 2 , the two planes are locked 
and the pinning length is larger than in a single plane system with the same disorder 
strength. This situation is reminiscent of the one that obtains in higher dimensions with a 
solid phase without dislocation at weak disorder and a more strongly pinned solid phase with 
dislocation at higher disorder. Such situation is expected on physical grounds, since forming 
dislocations allows the system to gain more energy from the random pinning potential. 

At non-zero temperature, the solid phase with dislocations is more easily melt than the 
solid without dislocations. Therefore, one expects that the pinning length decreases faster 
with temperature in the phase in which the vortices are pinned independently in the planes 
than in the phase in which the vortices are pinned by the antisymmetric component of 
disorder. Such behavior is indeed what we have obtained for T > 0. 



V. CONCLUSION 

In this paper, we have investigated two chains of strongly correlated bosons, coupled 
by boson interchain hopping and boson-boson repulsion, and applied our results both to 
quantum systems and to their classical analog of vortices in type II superconductors in the 
presence of columnar defects. 

Using bosonization and renormalization techniques, we have obtained the phase diagram 
in the absence of disorder. It contains four different phases: two charge density wave phases 
of different periodicity, a conventional superffuid phase, and a superfluid phase made of pairs 
of bosons. This last superfluid phase is an analog of the an XY phase found in spin one@ 
chains and spin ladderscj. For weak inter chain interactions one finds a transition between a 
superfluid and "crystalline" (charge density wave) phase. This transition occurs for consider- 
ably more repulsive interactions (K = 1/2) than in the case of a single chain, where it would 
occur for hard core bosons (K = 1). Interchain hopping thus stabilizes considerably the 
superfluid phase. Paradoxically, the long wavelength part of the interchain repulsion favors 
hopping between the chains and thus also helps in stabilizing a superconducting phase. 

We have then considered the effects of a random potential. For a single chain this leads 
to the localization of the bosons for moderately repulsive interactions (K = 3/2). In the 
case of the ladder the superfluid phase is much more resistant to localization. The transition 
occurs at more repulsive interactions (K = 3/4). In particular bosons with only contact 
interactions are now always superfluid. Two different localized phases exist: first a pinned 
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charge density wave where interchain hopping occurs, then for more repulsive interactions 
(K < 1/2), a pinned charge density wave where hopping is frozen and with "crystalline" 
order of the bosons. Interchain boson repulsion pushes the system towards localization. Us- 
ing renormalization we have also extracted the localization length as a function of disorder 
strength and interactions, and the temperature and frequency dependence of the conductiv- 
ity. For the ladder, the exponent of the superfluid correlation function takes at the transition 
the universal value rj = 1/3. The results for the compressibility and exponents obtained in 
our microscopic model are in agreement with the general scaling relations proposed for the 
superfluid-insulator transitionta. Moreover, the superfluid correlation exponent r\ for the 
ladder takes the same value than the one for a single chain@, suggesting the interesting 
possibility of a universal exponent of the superfluid-insulator transition in one dimension. 
Such a property would be much stronger than the simple inequality r\ < 1 suggested by the 
scaling. For an even number of chains, we have given a general argument, valid when all the 
mode but one are gapped, that this is indeed the case. This conjecture remains of course 
to be tested in much more details, both for the case of an odd number of chains and for 
situations where some modes could remain gapless. 

These results have a direct translation for the case of vortices in type II superconduc- 
tors. The bosonic ladder corresponds here to two planes of vortices with both vortex-vortex 
interactions and hopping of vortices between the planes. In the boson vortex mapping super- 
fluidity of the boson corresponds to a liquid state for the vortices, whereas the CDW phases 
correspond to crystalline order of the vortices. For the pure systems there is a transition 
between a high temperature liquid phase and a low temperature solid phase. The melting 
temperature T m is much smaller than the melting of a single plane due to the possibility 
of forming superkink between the planes. The transition is direct without an intermediate 
"supersolid" phase having both crystalline and superfluid properties, advocated to exist in 
this geometryB. This absence of supersolid phase could be due to our restriction to only 
two coupled planes instead of an infinite number. However since the existence of the su- 
persolid was only advocated based on perturbative arguments that would also give such a 
phase for the ladder, our study of the ladder shows that the possible existence of such a 
phase needs to be put on a firmer basis. In addition to these two phases, stable for weak 
interplane interaction, two other phases can exist: a vortex solid but where the vortices form 
superkinks between planes and a bound vortex pairs liquid. The random potential in the 
bosonic system is equivalent to columnar defects for the vortices. In presence of disorder two 
different transitions take place. At low temperature one is in a pinned solid state with no 
hopping between the planes. This corresponds to a Bose glass state for the vortices without 
dislocations between planes (nicknamed Bragg Bose glassful). Upon increasing the temper- 
ature a first transitions occurs at T m between the Bose glass and a phase where vortices 
hop between planes giving a (pinned) solid with half the periodicity of the low temperature 
one. Combining RG with Larkin-Ovchinnikov arguments allows to compute the critical cur- 
rents. At this transition the critical current is strongly decreased but remains finite. This 
decrease of the critical current is typical of columnar disorder and would not occur for pin- 
ning by twins. The liquid, with zero critical current is only recovered at higher temperature 
T = 3/2T m both for twins and columnar defects. This increase of the melting temperature 
is reminiscent of what is expected for columnar defects in higher dimensionality. Note that 
in our one dimensional system even an infinitesimal disorder would cause a finite shift of the 
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melting temperature. 

If disorder is increased, the planes decouple or equivalently dislocations appear between 
them. One recovers in this case the properties of individually pinned planes. This leads 
to an increase of the critical current similar to the one occurring for point like disorder 
when going from a Bragg glass state to a vortex glass stated. The transition occurs here 
between a dislocation free Bose glass (Bragg Bose glass) and a conventional Bose glass. Upon 
increasing the temperature there is then a single melting to a liquid phase but at a much 
higher temperature T' = 3T m . One thus expects for correlated disorder the same topology 
of phase diagram as for point like disorder, with two different glass phases but with the 
difference that increasing the correlated disorder raises the melting temperature. Whether 
the superkink phase close to the liquid for weak disorder exists in a three dimensional system 
is an interesting question that would deserve more investigations. 

The ladder geometry thus provide an interesting starting point to study competition 
between elasticity and disorder in vortex systems while retaining realistic feature. It would 
be interesting to investigate the effect of point disorder and the interplay of columnar and 
point disorder in that system. This has already been done for a single plane systemS and it 
would be interesting to know how interplane hopping modifies the response of the system to 
competing disorder. Also, systems with more than two coupled planes could provide further 
insight in the properties of coupled vortex planes and may be relevant to experiments. 
Finally, a mean field theory investigation of coupled planes would help in comparing the 
results with experiments on vortex pinning in BiSSCO. 
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APPENDIX A: BOSONIZATION OF BOSONS 



In this section the semi-phenomenological bosonization procedure introduced in Ref. [IT] 
to describe the low energy physics of one dimensional interacting bosons is reviewed briefly. 
The reader already familiar with this technique should read it quickly to get acquainted 



with the notations. In Ref. ^ ^ is argued that a one dimensional interacting boson system 
always has an associated bosonized Hamiltonian describing the low energy, long momentum 
phonon excitation of the system : 



H= ^ 



dx 
2^ 



7/ 

uK ^ Y[ f + -{d^) 2 



II and satisfy canonical commutation relations : [II(a;), 4>(x')} = %8{x 



x 



(Al) 



According to 



ref. the operators of the original bosonic theory can be expressed as a function of the 
fields entering the Hamiltonian ( |A1[ ) . The bosonized form of the boson annihilation operator 
in one dimension is : 
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oo 

ij) B {x) ~ e t0 ^ ( V e 2mWx)+npox) I (A2) 



where = 7r ^ rfx'n(x') and a is a short distance cutoff. 
The bosonized form for the density is : 



Ft A 00 

p B ( x ) = (p + ^) V e 2lm ^ + w) (A3) 

7T 



m=— cx) 



where po is the average boson density. Boson Green's functions and density density correla- 
tion functions are obtained by the standard methods of bosonization0'il@ giving: 



1 x ^ 

(rp e ine(x,T) e -m6(0,0)^ 1 1 X - 



+ U 2 T 2 

/ T ^ e t2n<f>(x,T) e -t2n<f>(0,0)\ _ f __}__\ 

\X 2 + U 2 T 2 J 



n 2 K 



(A4) 



The boson creation operator is also the order parameter for superfluidity and the 27rpo 
component of the density is the charge density wave order parameter. Thus, ( [A"4l) shows 
that in the one dimensional bosonic system there is no long range order (in agreement 
with the Mermin- Wagner theorenJ^E) , but only quasi long range order. In particular, the 
susceptibilities associated with these operators diverge in the long wavelength low frequency 
limit. When K increases, the system has a stronger tendency to superfluidity, whereas when 
K decreases it has a stronger tendency to charge density wave order. 

u, K still have to be determined as a function of the microscopic boson-boson interac- 
tions. To do so, one has to relate them to ground state parameters of the interacting system. 
In Galilean invariant systems , an exact relation between u, K holdfl and simplifies con- 
siderably the situation. This relation is : uK = — , where m is the mass of a boson. 
Knowing this, only one relation between u, K and a parameter of the ground state suffices 
to determine u and K. A convenient quantity is the compressibility x defined by: 

(A5) 

E being the ground state energy of the system. The compressibility can also be expressed 
as a function of u, K in the form: 

X = — 2 (A6) 

and permits a determination of u, K as a functions of the parameters of the microscopic 
Hamiltonian. x contains all interaction effects for Galilean invariant systems. On physical 
ground, the higher the boson-boson repulsion, the lower the compressibility. Therefore, K 
decreases when repulsion increases. According to ( |A4[ ) this means that repulsive interactions 
tend to favor charge density wave order versus superfluidity. In the case of bosons interacting 
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with a repulsive ^-function potential, the ground state energy and the compressibility can 
be obtained exactly as a function of the interaction^. This allows for the determination 
of u, K. In the case of a more general interaction potential, one has to rely on numerical 
techniques such as Quantum Monte Carlo simulations!!. 

Bosons on a lattice lack Galilean invariance, so that uK ^ Nevertheless, ( |A6| ) 

remains valid. To determine u, K one needs another relation linking_u, K to a parameter 
of the ground state. One convenient parameter is the charge stiffness^ defined for a closed 
ring of perimeter L by: 

L ( <PE(<p) 
2 ^ 

if being a flux threading the system. For bosons that do not carry electrical charge, it is 
nevertheless possible to define a "charge stiffness" by using the sensitivity of energy levels 
to boundary conditions since the only effect of the flux ip amounts to a change of boundary 
conditions. From the definition ( |A7| ), one obtains : 

V = uK (A8) 

In Galilean invariant systems, the charge stiffness is independent of the interactions. Com- 
bined with ( |A6| ) this equation allows to determine u, K as a function of the original param- 
eters of the Hamiltonian in the most general case. For instance, for hard core bosons on 
a lattice, using an exact mapping onto the XXZ spin 1/2 chain in a magnetic field^ one 
can obtain an exact expression of u, K using the Bethe Ansatz solution of the spin chain 
problem. The more general case of soft core boson needs to be treated numerically. 




APPENDIX B: SUPERFLUID-DISORDERED EXPONENT 

Let us consider a ladder with p legs. The density on chain i can be expressed in terms 
of the field 4>i as 

~ e 2l(tH (Bl) 

As for the two leg ladder it will be convenient to introduce the various modes (equivalent to 
the symmetric and antisymmetric ones) between the chains. If one now assumes that every 
mode except the totally symmetric one is gapped (the fact that the totally symmetric one 
remains always gapless is a consequence of the Galilean invariance), one can then proceed 
as for the two leg ladder. The totally symmetric mode has the following expression: 

S = ^I> (B2) 
And all the modes can be decomposed as: 

<pi = + non-symmetric modes (B3) 
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Let us assume that all of the non-symmetric modes are gapped with an average value for 
the corresponding 0's fields. The coupling to disorder should only retain the massless mode. 
It is thus necessary to form higher order terms in perturbation theory to cancel all the other 
modes. It is easy to check that the remaining coupling to disorder is of the form 

H = [ dxV{x)e l2p/ ^ x (B4) 



since it is necessary to go to p-th order to get rid of all the other modes. (fB4|) obviously 
reproduces the two leg result. Using the same RG argument the transition now occurs 
for K = 3/(2p), i.e. more and more repulsive interactions. Similarly the superfluid order 
parameter on chain i is given by 

V>< ~ e ld * (B5) 

Using again the decomposition in the above modes and assuming that all gaped modes have 
a an average value for the 9 field, one gets 

ipi ~ e^ ds C (B6) 

where C is a universal constant coming from the massive modes. In that case the correlation 
function at the transition decays as 

<V*(rM(0)>~ (^f^ (B7) 
giving the universal exponent r\ — 1/3 regardless of the number of legs. 



APPENDIX C: DERIVATION OF THE HAMILTON! AN OF THE COUPLED 
PLANES IN THE TIGHT-BINDING APPROXIMATION 

The Hamiltonian for non-interacting bosons in a plane with two potential wells is : 

H = -^(V; + Vl) + V(y - Y,) + V(y - Y 2 ) (CI) 
V(x) represents an isolated potential well. Consider the Hamiltonian h: 

h = ~^l + V(y) (C2) 

and call (p(y) its eigenfunctions. In the tight binding approximation the eigenfunctions are 
sought in the approximate form : 

y) = 4>{x){a<f(y - Y 1 ) + (3ip{y - Y 2 )) (C3) 

and the Hamiltonian is of the form h + V. The Hamiltonian is simplified into a Hamiltonian 
with a hopping term: t± = J dytp(y — Yi)V(y — Y 2 )f(y — Y 2 ) In the case : V(y) = —V 5(y), 
it is possible to obtain the exact eigenfunctions of ( |C2"| ) and calculate t± as : 
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t± = -V Kexp (-« \Y 1 -Y 2 |) (C4) 

with k = This allows a rewriting of the Hamiltonian in the second quantized represen- 
tation as: 



H = dx 



/• 



- ^ V ^i(*) + - ^ V ***(z) + ^(^(^^(x) + H. c.) 



(C5) 



The wavefunction (second-quantized) being: ^(x, y) = ^i(x)ip(y — Y±) + \p2(£)<£>(y ~~ ^2)- 
Adding interactions leads to the following expression for intrachain interaction: 

and for interchain interaction: 

u(x - x') = J dydy'U(x -x',y- y') \ V {y - Y x ) | 2 | ip{y' - Y 2 ) \ 2 (C7) 

There are also interactions which do not conserve the number of bosons in one chain However, 
these interactions merely correct the hopping term. In the case of the delta-function potential 
well, and a delta function repulsion between vortices in the form: 

U(x -x',y- y') = U 5(x - x')5(y - y') (C8) 

u, v have the following expressions : 

u(x - x') = U - 

v (x - x ') = [e" 4K|ri - y21 +2k\Y 1 -Y 2 \ e" 2 ^-^] (C9) 
8 l 

The tight binding approximation is well justified in the limit k \ Z\ — Z 2 |>> 1 only. Therefore 
in the case of delta functions interactions, interchain hopping dominates over interchain 
repulsion. This may change in the case of longer range interactions. 

APPENDIX D: DERIVATION OF THE COUPLING WITH DISORDER FOR 

THE TWO PLANE SYSTEM 

In the mapping of vortex planes onto coupled bosonic chains, introduction of columnar 
disorder in the two plane system is equivalent to introduction of a time independent random 
potential in the two chain system. The part of the Hamiltonian describing such disorder is: 

^disorder = J dxdyv(x, y)^(x, y)V(x, y) (Dl) 

with: 



v(x, y)v{x', y') = DS(x - x')S(y - y') (D2) 
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y being the coordinate perpendicular to the chains and \I/ being the full second quantized 
wave function. The chains being treated in the Tight Binding Approximation, one has: 

¥(ar, y) = ipx{x)ip{y - Y 1 ) + ik(x)<p(y - Y 2 ) (D3) 

Porting ( P3|) into ( |D1|) leads to the following coupling: 

^disorder = J dx (vi(x)lp{ljji(x) + V 2 (^2^2 (>)) + / dx5t ± (x) (^2 (x) + ififafa)) ( D4 ) 

With: 

M x ) = J dyv(x,y)tp 2 (y - Y 2 ) 
vi{x) = J dyv(x, y)(f 2 {y - Y x ) 
t±(x) = J dyv(x,y)ip 2 (y - Y 2 )<f> 2 (y - Y x ) (D5) 

These definitions lead to: 

V\(x)v\{x') = v 2 (x)v 2 (x') = D5(x — x) J dy(p A (y) 



vx(x)v 2 (x') = D6(x - x') J dy V 2 (y - Y 1 ) l p 2 (y - Y 2 ) 

6t ± (x)6t ± (x') = D5(x - a/)<£(y)<4(v) ( D6 ) 

If the potential that confines vortices in the plane is of the form —Vo(5(y — Y\) + S(y — Y 2 )), 
one can obtain an explicit form of the integrals in (p6|) . Introducing k = n ^ L , one has 
f(y) — \fi^^~ Ky and : 

f dyip 4 (y) = £ 

J — oo Z 

/OO fcf] 
dy V \y - YiV (V - Y 2 ) = «(1 + -)e~ 2 ™ 
-oo Z 

(D7) 

with a =| Zi — Z 2 | . The tight binding approximation being valid only for na ^> 1, V\ and 
v 2 are almost uncorrelated. Also, disorder induces a very weak random interchain hopping. 
For the tight binding approximation to be valid, the random interchain hopping has to be 
much smaller than t± . More precisely, the correlation length induced by the random hopping 
terms has to be much larger than the correlation length induced by t±. We know that : 
t± = — KVoexp ( — k | Z\ — Z 2 |). This leads to the following criterion: 

2 1 
1 \ i-Ka ( 1 



tj « bJ (D8) 

In terms of Vq the criterion reads : 

2(3-Kq) 

D « ^ "L (D9) 
k{1 + f )e< 4 -^) 

Since na ^> 1 when the tight binding approximation is correct, the condition on D is always 
satisfied as long as the tight binding approximation is correct. 
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APPENDIX E: RELATION BETWEEN ELASTIC CONSTANTS OF A ONE 
VORTEX PLANE SYSTEM AND THE PARAMETERS OF ITS EQUIVALENT 

BOSONIC THEORY 



We consider a plane containing vortices. Let p the vortex density, and c 44 and c\ the 
elastic constants of the vortex system. The elastic Hamiltonian of the vortex system is: 



H= 1 - 
2 



dxdz 



c AA (d z uf + ci(d x u) 



The Euclidean action of the equivalent bosonic system is: 



S 



dxdr 
2ttK 



u 



+ u(d x <f>f 



(El) 



(E2) 



Identifying the path integrals, and reminding that <fr(x) = 7ip u(x) one finds: 



C44 — T 



ci = T 



uK 
K 



(E3) 



Ci is the inverse compressibility both of the vortex system and the boson system. Also, for 
a translationally invariant vortex system: C44 = p§m Finally, K = 



APPENDIX F: ABSENCE OF A SUPERSOLID PHASE IN THE TWO PLANE 

VORTEX SYSTEM 

The discussion is done one the equivalent boson Hamiltonian. At the point K a = 1/2, a 
refermionization of H a can be performed!!!. One obtains: 



H a = [ ^v F [(vrn) 2 + (40) 2 1 + ^ I cos(2&) + -^r f cos(20 a ) (Fl) 

Where: a = V2<j) a and 6 a = y/29 a . 

Introducing the fermion fields f = e'( e °-+^) Qne obtains ^ e fermionized 
Hamiltoniar£3: 

H = -iv F J dx(¥ R d x ^ R - tf L d x V L ) + 2t ± J dx{¥ R ¥ L + m L ^ B ) + -jLJ dx(tf R V L + ^ R 

(F2) 

Such Hamiltonian can easily be diagonalized@. It is clear that expectation values such as: 
(^r) are zero on the ground state of this Hamiltonian due to the Mermin- Wagner theorem 
. This implies (sin( v / 20 a ) cos( v / 2^)) = 0. If one could have (9 2 a ), (0 2 ) 

,< 00 one would 

expect such quantity to be non-zero. This indicates that the presence of a supersolid phase 
in the two chain system is impossible. In a system of many coupled chains, however, the 
Mermin- Wagner does not apply and we cannot rule out a possible supersolid phase. 
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TABLES 



TABLE I. the 4 sectors of the two coupled bosonic chains model 





I 


II 


III 


IV 


K s 


< 1 


< 1/4 


> 1 


> 1/4 


K a 


< 1/2 


> 1/2 


< 1/2 


> 1/2 


a ,<f>a 


= 7=8 


(0a) = 


(4>a) = 


(0a) = 


phase 


CDW- 


CDW 47r P° 


BPSF 


SF+ 


Order Parameter 








e V2 
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FIGURES 




FIG. 1. The schematic phase diagram in the antisymmetric degrees of freedom as a function 
of K a and t±. The line separates the phase with (<^) < oo from the phase with (0%) < oo. For 
K a > 1 the system cannot sustain a phase with (<^) < oo. For 1/4 < K a < 1 a nonzero t± is 
needed in order to stabilize a phase with (9%) < oo. For K a < 1/4, the only possible phase is the 
phase with (4> a ) = 



CDW 4 



7ip f 
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FIG. 2. The phase diagram of the two chain bosonic system without disorder. The dashed line 
separates the domain of dominant CDW 47rpo from the domain of dominant SF+ fluctuations. The 
separation between dominant CDW- fluctuations and dominant BPSF fluctuations is at K s = 1. 
the dotted line corresponds to weakly coupled bosonic chains. 
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FIG. 3. The phase diagram of the two chain bosonic system in the presence of a random 
potential. The dotted line corresponds to weakly coupled bosonic chains. The hashed zones 
correspond to disordered phases in which bosons are localized. These disordered phases are pinned 
charge density waves. The stabilization of charge density waves by disorder permits to have a 
disordered 4-irpo charge density wave for the two weakly coupled chains. 
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FIG. 4. The two planes vortex system. The vortices are confined in the planes and can pin 
on columnar defects parallel to the magnetic field. In the situation represented on the figure, the 
temperature is low and the vortices are confined to the planes 
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(b) (d) 

FIG. 5. The behavior of vortices in the 2irpQ solid and in the 4irpo solid, (a) In the 2irpQ solid, 
the vortices are confined in the planes and the average distance between two vortices is 1/po- 

(b) Top view of the 2npo solid. 

(c) In the 47T/9o solid, vortices are hopping from one plane to the other, leading to a halving of the 
distance between two maximums of the vortex density. 

(d) Top view of the 47r/Jo solid. 
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Solid 



Liquid 







FIG. 6. The phase diagram of the pure two coupled planes vortex system. For T < T m one 
has the solid phase, with vortices confined in the planes, and for T > T m , one has the liquid phase 
with vortices hopping from plane to plane. 
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FIG. 7. The phase diagram of the weakly disordered two coupled planes vortex system. For 
T < T m , one has the pinned solid phase, for T m < T < Tl = §T m , one has the entangled pinned 
solid, in which vortex pinning is weaker and vortices can hop from plane to plane. For T > Tl, 
one goes into a vortex liquid phase in which vortices are completely unpinned. 
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FIG. 8. Dependence of critical current upon temperature for D = 0.1. For T > T m , entangl 
:nt appears and the critical current 
melt leading to a zero critical current. 



ment appears and the critical current is strongly reduced. For T > §T m , the entangled solid 
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FIG. 9. Phase diagram as a function of temperature and disorder. For large disorder, a Bose 
glass with dislocations is formed. At weak disorder, depending on temperature, there is a Bragg 
Bose glass (without dislocations), an entangled Bose Glass and a vortex liquid. The shaded squares 
correspond to regions where a multicritical point may lie. 
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